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ABSTRACT 

The small separation (<5^oi, 8vq2 and Sviz) between the oscillations with low degree 
I is dependent primarily on the sound speed profile within the stellar core, where 
nuclear evolution occurs. The detection of such oscillations for a star offers a very 
good opportunity to determine the stage of its nuclear evolution, and hence its age. In 
this context, we investigate the Sun and a Cen A and B. For a Cen A and B, each of 
the small separations <5^oi, <^02 and 5v\z gives a different age. Therefore, in our fitting 
process, we also employ the second difference, defined as v n ,2 — 2z/„ 1 4- v n ,o, which 
is 2<5^oi — 5vq2- In addition to this, we also use frequency ratio (v n ,0 / v n,2) ■ For the 
Sun, these expressions are equivalent and give an age of about 4.9-5.0 Gyr. For a Cen 
A and B, however, the small separation and the second difference give very different 
ages. This conflict may be solved by the detection of oscillation frequencies that can be 
measured much more precisely than the current frequencies. When we fit the models 
to the observations, we find (i) Z Q = 0.020, t = 3.50 Gyr and M B = 1.006M Q from 
the small separations Svqi, 5vq2 and 5v\3 of a Cen B; and (ii) a variety of solutions 
from the non-seismic constraints and S1/Q2 of a Cen A and B, in which the masses of a 
Cen A and B are slightly modified and the age of the system is about 5.2-5.3 Gyr. For 
Z = 0.025, the closest masses we find to the observed masses are Mb=0.922 M q and 
Ma=1.115 Mq. The differences between these masses and the corresponding observed 
masses are about 0.01 M Q . 

Key words: Sun: fundamental parameters- Sun: interior - stars: evolution stars: 
individual: a Cen - stars: interior 



1 INTRODUCTION 

Microscopic or macroscopic, solid or fluid of any kind, ev- 
ery object oscillates at frequencies determined by its struc- 
ture. In many cases, the oscillations sound the structure 
of the visible and invisible parts of the object in which 
they are trapped. Among such objects are stars, particu- 
larly solar-like stars, which oscillate in a variety of modes so 
that one can infer physical conditions deep in the nuclear 
core (Tassoul 1980) , which is not directly observable. In this 
context, Christensen-Dalsgaard (1988) developed a seismic 
Hertzsprung-Russel (HR) diagram in which one can deduce 
the size and evolutionary phase (age) of a star from its os- 
cillation frequencies. In this diagram, the horizontal axis is 
so-called the large separation between frequencies of consec- 
utive oscillations with order n (Af n j — v n ,i — v n -i,i) and the 
vertical axis is the small separation between frequencies of 
oscillations with different harmonics. As a result of nuclear 
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evolution, the position of a star in the seismic HR diagram 
changes with time. Roxburgh & Vorontsov (2003, hereafter 
RV2003) argued that the ratio of the small separation to the 
large separation is more sensitive to time than the accurate 
value of the small separation, and proposed an alternative 
expression for it. 

In many respects, and particularly because of its high 
quality data, the Sun is an excellent object to study stel- 
lar interiors for calibration of evolution code to be useds. In 
other words, our success in modelling interiors of stars of 
different kinds depends on how successful we are in mod- 
elling the solar interior. The solar models with chemical 
composition given by Grevesse & Sauval (1998) are in very 
good agreement with seismic inferences (Bahcall et al. 1995; 
Christensen-Dalsgaard et al. 1996; Gabriel & Carlier 1997; 
Yildiz 2001). However, the agreement disappears if the re- 
cent solar chemical composition given by Asplund, Grevesse 
& Sauval (2005, hereafter AGS2005) is used (see Basu et 
al. 2007 and references therein). Although uncertainties in 
the abundances of heavy elements are discussed by Pinson- 
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neault & Delahaye (2006) and found large enough to cover 
the old solar composition, an increase in diffusion coefficient 
(Guzik, Watson & Cox 2005) and opacity (Bahcall et al. 
2005) below the convective zone are suggested to restore the 
agreement between the models and helioseismic constraints. 

In our previous paper on a Cen (Yildiz 2007), we ob- 
tained 8.9 Gyr for the age of the system from the non-seismic 
constraints (NOS models in Yildiz 2007) and about 5.6-5.9 
Gyr from the seismic constraints (SIS Models in Yildiz 2007; 
see also Miglio & Montalban 2005). There may be obser- 
vational or theoretical reasons for this discrepancy. In this 
study, we investigate how well core structures of the Sun and 
a Cen A and B are represented by different forms of small 
separations and how precise their fundamental properties 
are. 

As main-sequence (MS) stars evolve their oscillation fre- 
quencies decrease because of the increase in their size. How- 
ever, the amount of decrease in the frequencies depends also 
on how the sound speed changes in the cavity in which the 
oscillations are trapped. As a result of nuclear evolution, the 
sound speed gradually decreases in the nuclear core, and 
therefore frequencies of oscillations with I = are much 
more reduced than those are with I — 1 and / = 2. Thus, 
the frequency ratio f n ,o/^n,2 is a function of time and has in 
principle diagnostic potential as the small separation for in- 
formation about the age of solar-like stars. We also use this 
expression, despite its lower sensitivity to time than 5vo2, in 
order to test if it gives the same age as the usual expressions 
(see Section 3). 

In Yildiz (2007), the small separation Do is computed 
from 8vo2 for a Cen A, and from <5^o2 and 5vi3 for a Cen B 
(Do = (5^02/6 + 5i/i 3 /10)/2), since 5v i = (v-no - (fn-1,1 + 
v n i)/2) values of models are quite different from the values 
found from the observed frequencies. In this study, using 
in addition the seismic data of a Cen A given by Bedding 
et al. (2004, BK2004), we consider if it is possible to fit 
model values of Svoi/2, 8V02/6 and Svis/W one by one to 
the corresponding values inferred from the seismic data of 
the Sun and a Cen A and B. 

The remainder of this paper is organized as follows. In 
Section 2, the basic properties of solar models with old and 
recent chemical composition are presented. The results are 
presented and discussed in Section 3. Finally, we make some 
concluding remarks in Section 4. 



2 PROPERTIES OF THE SOLAR MODELS 

The characteristics of the code used in the construction of 
models for stellar interiors are given in Yildiz (2007). Further 
details of the code are given in the references of that paper. 

Models of a Cen A and B required for our analysis 
are given in Yildiz (2007). We construct new solar models, 
because significant changes have occurred in our knowledge 
of the chemical composition of the Sun (AGS2005). 

The basic properties of the solar models with the re- 
cent solar composition given by AGS2005 are listed in Ta- 
ble 1. In the first row of this table the solar model with 
Zo = 0.016, which is in very good agreement with the ob- 
served chemical composition of the Sun, is presented: its 
surface helium and heavy element abundances are 0.244 and 
0.0124, respectively. However, the base radius of its convec- 
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Figure 1. The relative sound speed difference between the so- 
lar models and the Sun is plotted with respect to the relative 
radius. The first three models in the label arc constructed with 
the chemical composition given in AGS2005. The solid line with 
diamonds and dotted line show the solar models with Zo = 0.016 
and Zq = 0.017, respectively. The thick solid line represents the 
solar model with Zq = 0.016 and t© = 4.9 Gyr (see Section 3.4). 
The thick solid line with the filled circles is for the solar model 
with the chemical composition given by Grevesse & Sauval (1998), 
which is in good agreement with the helioseismic results. 



tive zone (0.733R©) is significantly greater than the value 
inferred from helioseismolgy (O.713±O.OO1R0 , Basu & Antia 
1997). In order to test the influence of initial value of heavy 
element abundance, another solar model is constructed with 
Zq = 0.017. This model is given in the second row of Table 1. 
The base radius of its convective zone (O.73OR0) is slightly 
smaller than that of the solar model with Zo = 0.016, but 
an agreement with the value inferred from helioseismolgy is 
not achieved. The age of these two models is taken as 4.6 
Gyr. For later usage, we also construct a solar model with 
Zo = 0.016 and age = 4.9 Gyr (see Section 3.4). 

In Fig. 1, the relative sound speed difference between 
these solar models and the Sun is plotted with respect to the 
relative radius. The relative sound speed difference between 
the solar model with Zo = 0.016 and the Sun is represented 
with diamonds. The largest difference, about 1.7 %, occurs 
in the region just below the base of the convective zone. 
The larger the initial heavy element abundance, the smaller 
the relative sound speed difference in the outer part of the 
radiative interior. For Zo = 0.017, for example, the largest 
difference reduces to 1.4 % (the dotted line). A very similar 
profile for the relative sound speed difference occurs for the 
model with Zo = 0.016 and age = 4.9 Gyr (solid line). For 
comparison, the model with Zo = 0.02 and mixture given 
by Grevesse & Sauval (1998) is also plotted in Fig. 1. It is 
in better agreement with the helioseismic results than the 
solar models with AGS2005. 
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Table 1. Properties of the solar models with the chemical composition given by AGS2005. The last line is for the observed values. The 
uncertainty in R c is 0.001 Rq (Basu & Antia 1997). The observed values of the surface Helium and metal abundances are taken from 
Basu & Antia (1995) and AGS2005, respectively. 
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Figure 2. The turning points of the observed modes in a Cen A 
and B computed from SIS models of Yildiz (2007). 



3 STRUCTURE OF NUCLEAR CORE AND 
OSCILLATION FREQUENCIES 

The acoustic oscillations are trapped in the stellar cavity. Di- 
mension of this cavity depends on degree (I) and frequency 
of oscillations. Whereas the cavity for the oscillations with 
I — extends down to the center, the oscillations with higher 
I have shallower cavities. The turning point (r t ) of an oscil- 
lation with non-zero I is approximately given as 



(1) 



Ct UJ 

where ui is the frequency of oscillation and c(r t ) is the sound 
speed at the turning point. The turning points of the oscilla- 
tions with observed frequencies of a Cen A and B are com- 
puted from SIS models and plotted with respect to n for the 
modes with I — 1 and I = 2 in Fig. 2. Whereas the modes of 
a Cen A (thin solid line with diamonds) and B (thin dotted 
line with xs) with I = 2 turns at nearly rt/R* — 0.1, the 
modes with I = 1 (thick solid line with boxes for a Cen A 
and dotted line with filled circles for a Cen B) sink deeper 
into the nuclear core (rt/R* ~ 0.05). 



Figure 3. The small separation Svo2 between the acoustic oscil- 
lations of the Sun (diamonds), the components of a Cen A (boxes 
for BC2002, filled circles for BK2004 and triangle for FC2006) and 
B (circles for Kjeldsen et al. and star for Carrier & Bourban 2003) 
as a function of n. 



Whereas the small separation is a sensitive function of 
physical conditions in central regions, where the nuclear evo- 
lution occurs, the large separation, Avo = v n ,o — ^n-i,o, 
is a measure of the mean density. The small separation 
between frequencies (Christensen-Dalsgaard 1988), 5i>02 = 
v n fi — Vn-i,2, is plotted as a function of n in Fig. 3, us- 
ing the observed frequencies of the Sun (thin solid line with 
diamonds; Chaplin et al. 1999), a Cen A (thick solid line 
with boxes; Bouchy & Carrier 2002, BC2002) and a Cen B 
(thin solid line with circles; Kjeldsen et al. 2005, hereafter 
KB2005). Its n dependence is very similar to a straight line 
for the Sun and a Cen A (BC2002). Although its average 
value is about 5.6 pHz for a Cen A, the difference between 
its maximum and minimum values is about 5 /jHz. In other 
words, the small separation has a strong dependence on n. 
Also shown in Fig. 3 are the small separations of a Cen A 
(triangle) and B (stars) computed from the observations of 
Fletcher et al. (2006, FC2006) and Carrier & Bourban (2003, 
CB2003), respectively. We note that the 8v l)2 from FC2006 
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Figure 4. Comparison of the small separation Suq2 computed 
from the old (dotted line with xs) and new (diamonds) acoustic 
oscillations of the Sun. 5vq2 of a Cen B (circles; KB2005) is also 
plotted. 

is significantly smaller than that of BC2002. The dotted line 
with filled circles in Fig. 3 is for 8vo2 from BK2004. Although 
the difference between its maximum and minimum values is 
about 10 /iHz, its mean value (6.3 /iHz) is close to the value 
found from BC2002 (5.6 /iHz). 

In Fig. 3, 5vo2 of the Sun (Chaplin et al. 1999) is much 
smoother than the corresponding values for a Cen A and B. 
However, for the less accurate solar data of the early stages 
of helioseismology this was not the case. In Fig. 4, Spq2 of 
the old (Grec et al. 1983; GF1983; dotted line with xs) and 
the recent (Chaplin et al. 1999) data of the Sun are plotted 
with respect to n. For small values of n, despite the scat- 
tering of the Grec et al. (1983) data, the two data sets are 
in agreement. However, for n > 25, there is a significant 
difference. Grec et al. (1983) state that the frequency reso- 
lution is 2 /iHz. The error in frequencies found by Chaplin 
et al. (1999) is about 0.05 /iHz for n ^ 22 and increases very 
rapidly for larger values of n: for n = 26, for example, it is 
about 0.3 /iHz. We also note in Fig. 4 the similarity between 
the old data of the Sun and the data of a Cen B. As our 
skill in detecting the seismic properties of a Cen A and B 
develops, as in the solar case, we may obtain much more 
smooth variation of Suq2 with respect to n. 

On the theory side, however, uncertainty arises from the 
fact that model frequencies are computed assuming an adi- 
abatic process for the oscillations (Christensen-Dalsgaard & 
Thompson 1997). The adiabatic approximation is valid for 
almost the entire interior, except for the near-surface re- 
gions. Therefore, the perturbation should also be applied to 
the energy equation, but the perturbation of energy equa- 
tion is highly uncertain, at least for the convective flux 
(Christensen-Dalsgaard et al. 1996). However, use of either 
the difference between or the ratio of frequencies for corn- 



Figure 5. Sound speed throughout the Sun and a Cen A and B 
with respect to the relative radius near the central regions. 

parison is a good way to minimize the effects of such a trou- 
blesome problem. 

3.1 Sound speed profiles of the Sun and a Cen A 
and B and the time variation of frequencies 

The variation of stellar oscillation frequencies with time is 
primarily a result of global expansion of the star like a 
heated ball. The secondary effect on the time dependence 
of frequencies arises from the fact that sound speed pro- 
file changes throughout the star. For the secondary effect, 
the most significant change occurs within the nuclear core. 
Therefore, this secondary effect is important for modes with 
low degree. Whereas variation of the large separation with 
time is a good measure of the first effect, the variation of 
the small separation with time is the result of the secondary 
effect. In Fig. 5, the sound speed profiles in the central re- 
gions of the Sun (thin dotted line, Basu et al. 1997), the 
solar model with Zo = 0.016 (stars), SIS models of a Cen A 
(thick dotted line) and B (thin solid line) are plotted with 
respect to the relative radius. For comparison, the sound 
speed (thick solid line) of a Cen A (SIS model) near the 
zero-age main-sequence (ZAMS) (t = 0.1 Gyr) is also plot- 
ted. As stated above, the most significant difference between 
the sound speed profiles occurs in the nuclear core: as hy- 
drogen is converted to helium, the mean molecular weight 
increases and consequently the sound speed at the center 
drops from 5.6 x 10 7 to 4.5 x 10 7 cm s _1 . 

The oscillations with I = 1 are influenced much more 
than the oscillations with I — 2 because of this variation 
in sound speed profile with time. It is well known that the 
turning point of oscillations depends mainly on the degree of 
oscillations. Whereas the turning points of oscillations with 
/ = 1 are about r t = 0.05i?*, the mean value for the turning 
points of oscillations with I = 2 is about r t = 0.10.R*. As 
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Figure 6. <5^oi2 as a function of n for a Cen A at various ages 
from the ZAMS to about the TAMS. For comparison, &vq\i of the 
solar model (Zq = 0.016 and age=4.6 Gy) and the SIS models of 
a Cen B are also plotted. 



Figure 7. The evolutionary tracks of a Cen A and B in the 
alternative seismic HR Diagram. For comparison, the observed 
position of the Sun and the position found from the solar model 
with Zq = 0.016 are also plotted. 



the nuclear evolution proceeds the sound speed decreases, 
and consequently the oscillations with I — 1 and I — 2 
sink deeper, according to equation (1). From the ZAMS 
to the present time, the decrease in sound speed at point 
r t = 0.05.R* with time for a Cen A, for example, is about 
0.6 x 10 7 cm s _1 (Fig. 5), at r t = 0.10.R*, however, the 
change with time is three times smaller (0.2 x 10 7 cm s _1 ) 
than this value. This means that the increase in the path of 
the oscillations with I = 1 is much more than that of the 
oscillations with I = 2. Therefore, the decrease in the fre- 
quencies of the oscillations with I = 1 as a result of nuclear 
evolution is much greater than that with 1 = 2. 

3.2 The second difference denned as 

As discussed above, oscillation frequencies of a star are a 
decreasing function of time, because of the increase of its 
size and sound speed profile over time as a result of nu- 
clear evolution. Christensen-Dalsgaard (1988) and RV2003 
use the frequencies of oscillations with different n and I for 
the expression for the small separation. We use also the sec- 
ond difference derived from the frequencies of mode with the 
same n: 

5^012 = (fn,0 - 2l/ n ,l + fn,2) (2) 

The second difference is the difference between 25i/oi and 
(5^02- According to the asymptotic relation, 5^012 is very 
small, just like 5vo2- In Fig. 6, <5^oi2 is plotted with respect 
to n at various evolutionary times of a Cen A, from the 
ZAMS to the terminal-age-MS (TAMS). Close to the ZAMS 
(thin solid line with diamonds), is negative and varies 

between —2.0 /kHz and —1.5 fiHz. As the model evolves, 



<Sfoi2 becomes positive at an age of nearly half of the MS 
life time of the star, and reaches an average value of about 
7.0 /iHz at age=7 Gyr. <5foi2 computed from the observed 
frequencies is also plotted in Fig. 6. From the comparison 
of the model and the observational 5voi2, contrary to the 
result of <5^o2 = v n fl — v n ,2, we deduce that a Cen A must 
be in the early phase of its MS evolution. 

The solar model with Zo = 0.016 and age=4.6 Gyr is 
also plotted in Fig. 6 (dotted line with circles). As for a Cen 
A, SP012 of the Sun is also independent of n. However, the 
model value of <5^oi2 for a Cen B is dependent on n near 
the ZAMS (solid line with xs) and the TAMS (thick solid 
line with stars) and has less n dependence between these 
two cases (dotted line with filled circles for a model with 
age=10 Gyr) . 

For the Sun, for which we have high quality of seismic 
data, average value of <5^oi2 derived from seismic data over 
n is —2.50 fiRz and it is —2.57 fiHz for the solar model with 
Zo = 0.016 and age=4.6 Gyr (dotted line with circles in 
Fig 6); that is, the observed value is slightly greater than 
the model value (the same result is obtained from Svo2 = 
v n ,o — v n ,2'i see Section 3.4). 

In Fig. 7, average value of <5foi2 over n is plotted against 
Auo during the MS evolution (SIS model) of a Cen A (solid 
line with diamonds) and B (dotted line with xs). The ZAMS 
is represented by the symbol closest to the legend '0.1 Gyr'. 
The next closest symbol corresponds to t = 1 Gyr. After 
that symbol, time increment between two symbols without 
legend is 1 Gyr. In the early phase of MS evolution, 8f i2 of 
a Cen B remains almost constant, whereas <5foi2 of a Cen 
A increases very rapidly. In contrast to 8vo2, again, Sf i2 of 
the massive star (a Cen A) is much more sensitive function 
of time than that of the low-mass star (a Cen B). 
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Figure 8. The ratio v n ,o/ u n,2 from SIS model of a Cen A with 
respect to n. For comparison, the observed ratio is also plotted 
(filled circles). 



3.3 Ratio of frequencies 

Although the ratio v n ,o/vn,2 is not a more sensitive function 
of time than 5vo2, we want to test whether it gives the same 
age as 5fo2- In Fig. 8, ratio of z/ nj o to ^ n ,2 for a Cen A is 
plotted as a function of n at 1 Gyr (thin solid line with 
diamonds), 3 Gyr (thin dotted line), 5 Gyr (thick dotted 
line) and 7 Gyr (thick solid line with boxes). The decrease 
in the ratio with time is noticeable. For comparison, the 
observed ratio is also plotted (solid line with filled circles). 
The observed ratio is in agreement with the curve of 5 Gyr, 
except for the ratio for n = 16. 



3.4 Age and metallicity of the Sun from the 
seismic constraints 

The observed value of 8vo2 for the Sun is found to be 9.84 ± 
0.1 fjRz from the seismic data of the Bison group. The value 
of &V02 for a solar model with Zo = 0.016 is, however, 10.11 
/xHz; that is, it is slightly greater than the observed value. 
For the solar model with Z = 0.017, 6v 2 = 10.08 //Hz. The 
value of 5vo2 for a solar model with Zo = 0.016 and age=4.9 
Gyr, 9.86 fiRz, is in very good agreement with the observed 
value. From the ratio of frequencies (t'n, 0/^11,2), we confirm 
that the same age (4.9 Gyr) is more suitable than 4.6 Gyr. 
Similarly, from the comparison of the observed value of <5^oi2 
(2.50 piHz) with that of the solar model with Zo = 0.016 and 
age=4.6 Gyr (2.57 //Hz), we find that the solar model with 
4.9 Gyr is in better agreement with the helioseismic data. 



3.5 Are values of 5uoi/2, 8P02/6 and 81/13/W 
compatible with each other 

The small separation Do can be derived from Svoi, 5vo2 and 
81/13, if we have enough observational data: 

D = (Doi+D 2+D 13 )/3 = (<5i/ O i/2+<^o2/6+M 3 /10)/3.(3) 

As emphasized above, in Yddiz (2007), Do is computed from 
SV02 for a Cen A and from 5^02 and 81/13 for a Cen B. Its val- 
ues inferred from the observed seismic frequencies (D 2, o bs 
and Di3,obs) and models with an age of about 5.7-5.9 Gyr 
are in good agreement. The values of Doi, bs = 5voi/2 of 
both a Cen A and B are not used in the analysis in Yddiz 
(2007), because they are not in agreement with the D01 val- 
ues derived from SIS models and are not compatible with 
Do2,obs and Di 3i0 bs (see Table 2). In fact, the disagreement 
between the values of <5^oi2 computed from the models and 
the values inferred from the seismic data (see section 3.2) 
arises from this inharmoniousness. Such an inharmonious- 
ness does not exist for the Sun, and we may hope to remove 
it with observational seismic data for a Cen A and B that 
are much more precise than the present data. 

We first consider the solar data and try to fit D01, 
D02 and D13 of the solar models to D O i, o bs0, D 2,obs© and 
Di3 jO bs0. In Table 2, D01, D02 and D13 of the solar models 
and models of a Cen A and B are listed. In the first three 
rows the solar models with an age of 4.60 Gyr and with dif- 
ferent Zo values are given. D01, D02 and D13 of these models 
are greater than D O i, o bs0, D O 2,ob S and Di3 i0 bs© (given in 
the eighth row), respectively. For comparison, the values of 
Do2,obs0 and Di3 jO bs0 derived from the recent seismic solar 
data of Basu et al. (2007) and their uncertainties are given 
in ninth and tenth rows, respectively. In the fourth row, the 
solar model with an age of 4.9 Gyr and Zo = 0.016 is given. 
Its D02 is the same as Do2,obs, and D01 and D13 are very 
close to D O i,obs0 and Di3 jObs 0, respectively. In order to find 
the optimum values of Zo and an age for the Sun we write 
three equations for each D u i , for instance, 



D, 



'01 ,obs0 



APoi ^ AD01 
Mn.ref H — A ^ 0Z0 H r^~ dt ' 



AZ ' At (4) 

and solve them for Zo and age (t) . We compute the required 
derivatives such as given in equation (4) from the solar mod- 
els discussed above and find the solution: Z = 0.017 and 
t = 4.93 Gyr. The solar model with these values is given in 
the sixth row of Table 2. We use this model as the refer- 
ence model and obtain another solution: Zo = 0.0165 and 
t — 4.98 Gyr. The solar model with the new results is pre- 
sented in the seventh row. D H i of these solar models with 
about Zo = 0.0165 - 0.0170 and age about 4.9 - 5.0 Gyr 
are very similar values and are compatible with the results 
inferred from helioseismolgy. 

In this context, we also consider a Cen A and B. Their 
models from Yddiz (2007) are given in the first three rows 
of their corresponding part in Table 2. For our analysis of 
a Cen A and B, we write again three equations similar to 
equation (4) for D tl i for each star but with three unknowns. 
In addition to Zo and t , we assume that the mass of the 
stars are also unknown. We obtain separate solutions for a 
Cen A and B. For a Cen B, Z = 0.020, t = 3.50 Gyr and 
M B = 1.006M Q . For a Cen A, we obtain t = 3.50 Gyr and 
Ma — 1.19Mq. However, value of Zo for the solution for a 
Cen A is not a reasonable value and therefore we use the 
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Table 2. Dqi, D02 and D13 of the solar models and models of a Cen A and B. 



Star 


t(Gyr) 


Z 


X 


Doi(mHz) 


D 02 (mHz) 


Di 3 (/iHz) 


Model 


Sun 


4.60 


0.0160 


0.70975 


1.879 


1.686 


1.739 






4.60 


0.0170 


0.70331 


1.870 


1.679 


1.734 






4.60 


0.0190 


0.69080 


1.858 


1.659 


1.713 






4.90 


0.0160 


0.71280 


1.825 


1.644 


1.701 






4.90 


0.0200 


0.68806 


1.825 


1.613 


1.676 






4.93 


0.0170 


0.70665 


1.837 


1.632 


1.695 






4.98 


0.0165 


71039 


1.835 


1.628 


1.691 




Obs C99 








1.813 


1.640 


1.679 




Obs B07 


_ 


_ 


_ 




1.653 


1.683 














± 0.005 


± 0.007 




a Cen A 


5.70 


0.0322 


0.669 


. 


0.93 


0.76 


SIS 




8.88 


0.0328 


0.703 


2.53 


0.49 


0.37 


NOSIS 




5.90 


0.0230 


0.713 


2.24 


0.94 


0.82 


SIS23 




3.50 


0.0200 


0.756 


3.01 


1.44 


1.11 


SIS2 




3.50 


0.0200 


0.732 


2.73 


1.41 


1.14 


SIS2p 


Obs BK2004 








1.41 


1.05 


1.10 




Obs BC2002 








0.98 


0.99 






Uncertainty 








0.12 








a Cen B 


5.70 


0.0322 


0.669 


0.95 


1.53 


1.86 


SIS 




8.88 


0.0328 


0.703 


-0.05 


0.92 


1.44 


NOSIS 




5.90 


0.0230 


0.713 


0.98 


1.49 


1.71 


SIS23 




3.50 


0.0200 


0.756 


1.83 


1.82 


1.81 


SIS2 




3.50 


0.0200 


0.732 


1.60 


1.82 


1.88 


SIS2p 


Obs KB2005 








2.26 


1.69 


1.67 




Uncertainty 








0.26 


0.06 


0.07 






same value of Zo as obtained for a Cen B. Models of a Cen 
A and B (SIS2) with these values are given in the fourth 
row of their corresponding parts in Table 2. Although these 
values are the optimum values in order to fit model values 
of D u i to Dn'.obs, there are still significant differences be- 
tween D01 and £>oi,o&s of each star. We construct similar 
models for a Cen A and B with the same values but the 
masses (Ma = 1.15M0 and Mb = 0.97M Q ) are average of 
the masses we find and the masses given by Pourbaix et al. 
(2002). These models (SIS2p) are presented in the fifth row 
of the corresponding parts of a Cen A and B in Table 2. D01 
values of the SIS2 and SIS2p models are significantly differ- 
ent. However, there is no simultaneous agreement between 
the model and the observed values of D H i for a Cen A and 
B as in the case of the Sun. 



4 MASSES AND AGE OF a CEN A AND B 

An agreement between the ages of a Cen A and B derived 
from their seismic and non-seismic constraints can be ob- 
tained by modification of their masses. To do this, we write 
down six equations (for luminosity, radius and D02 = Su 02 /6 
of the each component) similar to equation (4) and solve 
them to find Ma, Mb, qa, Qb, X and t for fixed Z. These 
equations are not independent and therefore we get a set of 
solutions. Some of these solutions are given in Table 3, for 
Z = 0.025 and Z — 0.020. For comparison, the values found 
from the observed data are also given (for a Cen A, Bazot 
et al. (2007) gives D = 1.15 ± 0.07/iHz). Different combina- 
tions of Ma and Mb are possible. In Fig. 9, Mb is plotted 
against Ma- The closest masses to the observed masses are 



0.9 1 1 1 1 1 1 1 

1.08 1.09 1.10 1.11 1.12 1.13 1.14 

M A /M 3 



Figure 9. By modification of masses of a Cen A and B, full 
agreement between the observational and theoretical results is 
achieved. However, there is no unique solution. Various combi- 
nations of Ma and Mb are possible. The mass of Mb required 
for the full agreement (circles) is plotted with respect to Ma for 
Zo = 0.025. The solutions for Zq = 0.020 are represented by filled 
circles. The observational value is also shown (boxes). 
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Table 3. Models with slightly modified masses, that are in full agreement with observations. 



Star 


M/M 


t(Gyr) 


Z 


a 


X 


L/L 


i?/R 


T cff (K) 




Avo(fiHz) 


A 


1.097 


5.15 


0.025 


1.588 


u . uoyuu 


1.544 


1.224 


5824 


1.072 


106.0 


B 


0.907 


5.15 


0.025 


1.550 


0.68950 


0.507 


0.864 


5243 


1.679 


162.3 


A 


1.113 


5.30 


0.025 


1.608 


<J. f U -LOO 


1.542 


1.228 


5809 


1.049 


106.1 


B 


0.921 


5.30 


0.025 


1.644 


0.70133 


0.508 


0.863 


5252 


1.670 


163.8 


A 


1.115 


5.24 


0.025 


1.611 


70256 


1.536 


1.225 


5812 


1.039 


106.6 


B 


0.922 


5.24 


0.025 


1.648 


0.70256 


0.505 


0.862 


5248 


1.679 


164.22 


A 


1.130 


5.25 


0.025 


1.634 


0.71180 


1.539 


1.224 


5816 


1.067 


107.4 


B 


0.934 


5.25 


0.025 


1.694 


0.71180 


0.507 


0.863 


5248 


1.687 


164.7 


A 


1.096 


5.35 


0.020 


1.575 


0.71830 


1.538 


1.223 


5855 


1.086 


106.0 


B 


0.910 


5.35 


0.020 


1.635 


0.71830 


0.508 


0.863 


5250 


1.681 


162.8 


A 


1.107 


5.37 


0.020 


1.582 


0.72504 


1.552 


1.228 


5819 


1.075 


105.7 


B 


0.918 


5.37 


0.020 


1.637 


0.72504 


0.507 


0.866 


5239 


1.680 


162.6 


A (obs.) 


1.105 










1.544 


1.224 


5824 


1.05 


105.6 


B (obs.) 


0.934 










0.507 


0.863 


5250 


1.69 


161.5 



Table 4. Adiabatic oscillation frequencies of the best fitted models for a Ccn A and B with masses 1.115 Mq and 0.922 M , respectively. 



a Cen A 



a Cen B 



n 






Vn2 








Vn2 


^n3 


10 


1219 


56 


12%7 


70 


1320 


00 


1108 


49 


1844 


75 


1924 


54 


1995 


22 


1676 


04 


11 


1327 


68 


1375 


06 


1426 


31 


1219 


56 


2010 


85 


2089 


51 


2160 


28 


1844 


75 


12 


1433 


86 


1480 


79 


1532 


34 


1327 


68 


2175 


42 


2253 


87 


2324 


48 


2010 


85 


13 


1539 


96 


1587 


38 


1639 


11 


1433 


86 


2339 


09 


2418 


09 


2489 


06 


2175 


42 


14 


1646 


69 


1694 


54 


1746 


41 


1539 


96 


2503 


03 


2581 


58 


2652 


87 


2339 


09 


15 


1753 


84 


1801 


49 


1852 


86 


1646 


69 


2666 


40 


2744 


84 


2815 


82 


2503 


03 


16 


1860 


08 


1907 


62 


1958 


69 


1753 


84 


2828 


74 


2907 


42 


2978 


97 


2666 


40 


17 


1965 


66 


2013 


16 


2064 


48 


1860 


08 


2991 


41 


3070 


09 


3142 


26 


2828 


74 


18 


2071 


30 


2119 


40 


2171 


02 


1965 


66 


3154 


33 


3233 


53 


3305 


99 


2991 


41 


19 


2177 


64 


2226 


17 


2278 


23 


2071 


30 


3317 


60 


3397 


04 


3470 


18 


3154 


33 


20 


2284 


65 


2333 


49 


2385 


51 


2177 


64 


3481 


40 


3560 


66 


3634 


05 


3317 


60 


21 


2391 


67 


2440 


74 


2492 


85 


2284 


65 


3644 


93 


3724 


53 


3798 


19 


3481 


40 


22 


2498 


74 


2547 


79 


2599 


94 


2391 


67 


3808 


69 


3888 


42 


3962 


70 


3644 


93 


23 


2605 


60 


2654 


97 


2707 


21 


2498 


74 


3972 


89 


4052 


78 


4127 


39 


3808 


69 


24 


2712 


60 


2762 


14 


2814 


69 


2605 


60 


4137 


29 


4217 


50 


4292 


41 


3972 


89 


25 


2819 


87 


2869 


70 


2922 


25 


2712 


60 


4301 


98 


4382 


13 


4457 


40 


4137 


29 


26 


2927 


18 


2977 


24 


3029 


99 


2819 


87 


4466 


71 


4546 


95 


4622 


34 


4301 


98 


27 


3034 


68 


3084 


85 


3137 


60 


2927 


18 


4631 


36 


4711 


65 


4787 


36 


4466 


71 


28 


3142 


05 


3192 


45 


3245 


21 


3034 


68 


4796 


13 


4876 


35 


4952 


35 


4631 


36 


29 


3249 


41 


3299 


82 


3352 


64 


3142 


05 


4960 


90 


5041 


26 


5117 


37 


4796 


13 


30 


3356 


62 


3407 


17 


3459 


85 


3249 


41 


5125 


69 


5206 


10 


5282 


43 


4960 


90 



M B =0.922 M and M A =1.115 M , for Z = 0.025. Whereas 
Mb is less than the observed value by 0.012 M , Ma is 
greater than the observed value by 0.010 M . We notice 
two common properties of these models: 

(i) thc age of the system is about 5.2-5.3 Gyr, 

(ii) for each solution, except the one given in the first and 
the second rows, the convective parameter (fifth column in 
Table 3) of a Cen A is smaller than that of a Cen B. This 
result is in good agreement with the result given in fig. 5 of 
Yildiz (2007). 

The large separations of the models of a Cen A and 
B whose radii are fitted to the observed radius are greater 
than the observed large separations. In order to fit the large 
separations of the models to the values inferred from astero- 



seismology, the radii of a Cen A and B should be J?a=1-230 
R and i? B =0.875 R . 

The adiabatic oscillation frequencies of the best fitted 
models for a Cen A and B with masses 1.115 M and 0.922 
M , respectively, are given in Table 4. The fundamental 
properties of these models are presented in the sixth and 
the seventh rows of Table 3. 



5 CONCLUSIONS 

Using solar models and models of a Cen A and B, the seismic 
properties of these stars are investigated in detail. Because 
of inconsistent results on the age of the a Cen system from 
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the classical and seismic constraints, we test how well the 
small separation indicates age. In this context, we consider 
two expressions in place of the customarily used form of the 
small separation between the frequencies (8vq2 or Do) : 

(i) the second difference, 5foi2 = (yn,o ~ 2v n ,i + ^11,2); and 

(ii) the frequency ratio (u ni o/vn,a) ■ 

For the Sun, these three expressions give consistent results: 
the solar models with age of 4.9-5.0 Gyr and with Zo « 
0.0165 are in better agreement with the helioseismic data 
than the solar model with Zq — 0.0165 and age 4.6 Gyr. For 
a Cen A and B the situation is complicated: although Svo2 
gives about 5.7-5.9 Gyr for the age of the system (nearly the 
TAMS of a Cen A), we deduce from the comparison of the 
observed and model values of 5^oi2 that a Cen A and B are 
not so much evolved. 

The expression we use for the small separation, <5^oi2, 
has several important advantages over the customarily used 
expression Svo2' 

(i) the resolution of the seismic HR diagram with <5i/oi2 for 
the massive solar-like stars (e.g., a Cen A) is higher than 
that with Svo2', 

(ii) the evolutionary track of a star in the new seismic HR 
diagram is in agreement with its evolutionary track in the 
classical HR diagram (luminosity-effective temperature); 
and 

(iii) 5foi2 is nearly independent of n. 

Si/012 deserves more detailed investigation in order to 
judge whether it has more specific advantages (Christensen- 
Dalsgaard, private communication). 

From the three equations for Dui — Suui / [Al + 6) for a 
Cen A and B, we obtain solutions for initial heavy element 
abundance, age and mass of each star: for a Cen B, Zo = 
0.020, t = 3.50 Gyr and Mb = 1.006M©; for a Cen A, 
t = 3.50 Gyr and M A = 1.19M . Z value found for a Cen 
A is not a reasonable value. Therefore, we adopt Zq — 0.020. 
The ratio of masses (Ma /Mb) is the same as the mass ratio 
derived from the observation (Pourbaix et al. 2002). The 
heavy element abundance at the surface of model of a Cen 
B with diffusion (SIS2) is Z s = 0.0179. This means that 
the average overabundance relative to solar (Z s = 0.0122, 
AGS2005) is 0.17 dex. 

Using the non-seismic constraints and 8vo2 of a Cen 
A and B we find from the solution of six equations that 
the age of the system is about 5.2-5.3 Gyr. We also confirm 
that the mixing length parameter of a Cen A is smaller 
than that of a Cen B. This is consistent with the result 
given in fig. 5 of Yildiz (2007) , which shows that the mixing 
length parameter of a Cen A is a decreasing function of 
time. Another interesting result we find is that the ages of 
the Sun and the a Cen system is very close to each other, 
about 5 Gyr. 

The models of a Cen A and B, that have the same lu- 
minosity and radius as the observed values, have large sep- 
arations greater than the observed values. If 7?a=1-230 Rq 
and i?B=0.875 Rq, the large separations of the models are 
in good agreement with the values inferred from asteroseis- 
mology. 

In comparing the old (Grec et al. 1983) and the recent 
(Chaplin et al. 1999) seismic solar data and find that the 
level of accuracy of the seismic data of a Cen A and B is 
very similar to that of the old solar data. We hope that the 



situation will be much better in the near future than the 
present. 
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